Theoretical approach to microwave radiation-induced zero-resistance states in 2D 

electron systems 
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We present a theoretical model in which the existence of radiation-induced zero-resistance states 
is analyzed. An exact solution for the harmonic oscillator wave function in the presence of radiation, 
and a perturbation treatment for elastic scattering due to randomly distributed charged impurities, 
form the foundations of our model. Following this model most experimental results are reproduced, 
including the formation of resistivity oscillations, their dependence on the intensity and frequency 
of the radiation, temperature effects, and the locations of the resistivity minima. The existence of 
zero-resistance states is thus explained in terms of the interplay of the electron MW-driven orbit 
dynamics and the Pauli exclusion principle. 
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PACS numbers: 

In the last two decades, especially since the discov- 
ery of the Quantum Hall Effect, a lot of progress has 
been made in the study of two-dimensional electron sys- 
tems (2DES), and very important and unusual properties 
have been discovered when these systems are subjected 
to external AC and DC fields. In the last two years two 
experimental groups^ have announced the existence of 
vanishing resistance in 2DES, i.e. zero resistance states 
(ZRS), when these systems are under the influence of a 
moderate magnetic field (B) and microwave (MW) radi- 
ation simultaneously. In the same kind of experiments 
large resistivity oscillations have been observediiSiii. 
The discovery of this novel effect has led to a great 
deal of theoretical activity, and among the most inter- 
esting contributions we can summarize various different 
approaches. Some, like Girvin et al;^ argue that this 
striking effect has to do with photon-assisted scattering 
from impurities or disorder, or alternatively arises from 
acoustic phonon scattering. OthergSiSi relate the ZRS 
with a new structure of the density of states of the sys- 
tem in the presence of light. According to Andreev's 
approach^", the key is the existence of an inhomogeneous 
current flowing through the sample due to the presence 
of a domain structure in it. To date there is no consen- 
sus about the true origin. All the theories above have 
in common that they predict negative resistivity, while 
this has not been experimentally confirmed. Only Wil- 
lett et alii, have recently observed negative conductivity 
for certain configurations of contacts. 
In this Letter we develop a semi-classical model that is 
based on the exact solution of the electronic wave func- 
tion in the presence of a static B interacting with MW- 
radiation, i.e. a quantum forced harmonic oscillator, and 
a perturbation treatment for elastic scattering from ran- 
domly distributed charged impurities. We explain and 
reproduce most experimental features, and clarify the 
physical origin of ZRS. For large MW field amplitudes, 
final states for electrons, which semi-classically describe 
orbits whose center positions oscillate due to the MW 
field, will be occupied. The MW field thus blocks the 



electron movement between orbits, and the longitudinal 
conductivity and resistivity pxx will be zero (see Fig.l). 
We explain the dependence of p^x on temperature (T) by 
means of the electron interaction with acoustic phonons 
which acts as a damping factor for the forced quantum 
oscillators, and we discuss the ZRS dependence on B. 
We observe them at w/wc = j-l-1/4 for the experimental 
parameters of ref . [1] , being w the MW frequency, Wc the 
cyclotron frequency and j an integer. In our model we do 
not consider the induced electrostatic potentials by the 
charge distribution within the sample, nor the dynamical 
electronic response induced by the AC-potential»i2iii 
We first obtain an exact expression of the electronic wave 
vector for a 2DES in a perpendicular B, a DC elec- 
tric field and MW radiation which is considered semi- 
classically. The total hamiltonian H can be written as: 
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X is the center of the orbit for the electron spiral motion: 



X = ^ 

eB 



ra'w^ ' -^0 the intensity for the MW field 
and Edc is the DC electric field in the x direction. Hi 
is the hamiltonian corresponding to a forced harmonic 
oscillator whose orbit is centered at X. Hi can be solved 
exactljii^ii^, and using this result allows an exact solution 
for the electronic wave function of H to be obtained: 
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where 7 is a phenomenologically-introduced damping fac- 
tor for the electronic interaction with acoustic phonons. 
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(j)n is the solution for the Schrodinger equation of the 
unforced quantum harmonic oscillator and Xci(t) is the 
classical solution of a forced harmonic oscillatosi^, Xd = 
, cos wt. L is the classical lagrangian, and 

Jm are Bessel functions. Apart from phase factors, the 
wave function for H is the same as the standard har- 
monic oscillator where the center is displaced by Xdit). 
Now we introduce the impurity scattering suffered by the 
electrons in our modeli&. If the scattering is weak we 
can apply time dependent first order perturbation theory, 
starting from H as an exact hamiltonian and '^i{x,t) as 
the wave- vector basis. The aim is to calculate the tran- 
sition rate from an initial state ^'„(a;,t), to a final state 
*m(a;,t): 
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where Vs is the scattering potential for charged 
impurities^"^: Vs = Eg 2Se(g+9.) ' ^ being the 

surface of the sample, e the dielectric constant and qs 
is the Thomas-Fermi screening constantii. After some 
lengthy algebra we arrive at the following expression for 
the transition rate: 
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With 



the experimental parameters we take, the arguments of 
the Bessel functions are very small (~ 10"^), and only 
Jo terms need to be considered. F is the Landau level 
broadening, rij is the impurity density and R is the mag- 
netic characteristic length R^ = L"i~"'2 g^j-g 
associated Laguerre polynomials, ni = max{n, m) and 
77.2 = min{n, m). 

Without radiation, an electron in an initial state cor- 
responding to an orbit center position scatters and 
jumps to a final state with orbit center X^^^, changing 
its average coordinate in the static electric field direction 
by = — X° = qcosOR^ (polar coordinates, q 

and 6, have been used). In the presence of MW radia- 
tion, the electronic orbit center coordinates change and 
are given according to our model by X^"*^ = + Xci{t). 
This means that due to the MW field all the electronic 
orbit centers in the sample oscillate back and forth in 
the X direction through Xd. We have to consider two 
important factors: first, when an electron suffers a scat- 
tering process with a probability given by Wm,n, it takes 
a time r 



for that electron to jump from an 
orbit center to another. Secondly, in the jump, as in 
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FIG. 1: Schematic diagrams of electronic transport without 
and with MW. In Fig.(l.a) no MW field is present, and due 
to scattering electrons jump between fixed-position orbits. 
When the MW field is on, the orbits are not fixed but os- 
cillate at w. In Fig.(l.b) the orbits move backwards during 
the jump, and on average electrons advance further than in 
the no MW case. In Fig.(l.c) the orbits are moving forwards, 
and on average electrons advance less than in the no MW 
case. In Fig. (l.d) the orbits are moving forwards but their 
amplitudes are larger than the electronic jump, and the elec- 
tron movement between orbits cannot take place because the 
final state is occupied. This situation corresponds to ZRS. 



any other scattering event, the electron loses memory 
and phase with reference to the previous situation, and 
thus when it arrives at the final state the oscillation con- 
dition is going to be different from the starting point. 
If we consider that the oscillation is at its mid-point 
when the electron jumps from the initial state, and that 
it takes a time r to get to the final one, then we can 
write for the average coordinate change in the x direc- 



tion: AX 



MW 



AX'' 
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In Fig. (1) we present schematic diagrams for the dif- 
ferent situations. In Fig. (l.a) no MW field is present 
and electrons jump between fixed orbits, and on average 
an electron advances a distance AA"". When MW field 
is on, the orbits are not fixed, and instead move back 
and forth through Xd. Three cases can be distinguished. 
In Fig. (l.b) the orbits are moving backwards during 
the electron jump, and on average, due to scattering pro- 
cesses, the electron advances a larger distance than in 
the no MW case, AX^^ > AA°. This corresponds to 
an increasing conductivity. In Fig. (I.c) the orbits are 
moving forwards and the electron advances a shorter dis- 
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tance, AX*^*^ < AX°. This corresponds to a decrease 
in the conductivity with respect to the case without MW. 
If we increase the MW intensity, we will eventually reach 
the situation depicted in Fig. (l.d) where orbits are mov- 
ing forwards, but their amplitude is larger than the elec- 
tronic jump. In that case the electronic jump is blocked 
by the Pauli exclusion principle because the final state is 
occupied. This is the physical origin of ZRS. At different 
ranges of parameters, i.e. larger MW power or smaller 
w, additional terms corresponding to Bessel functions of 
order higher than zero would contribute to Eq. 4, which 
may eventually produce negative conductivityii^i^. 
If the average value AX*^^ is different from zero over 
all the scattering processes, the electron possesses an av- 
erage drift velocity Vn^m in the x direction^^. This drift 
velocity can be calculated readily by introducing the term 
^j^MW jj-^^Q ^Yie integrand of the transition rate, and fi- 
nally the longitudinal conductivity be written as; 
cTxx = p{En)vnjm[f{En) ^ f{Em)]dE„. Gathering 
all the terms, we finally obtain the expression: 
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where the density of states p{En) has been simulated 
by a Lorentzian function, being A — 

To obtain px 
where Uxy — 



; we use the relation p^x — 
^ and ctt 



All our results have 



cui'a ^ XX ^xy- 

been based on experimental parameters corresponding to 
the experiments of Man:^ et al. In Fig. 2 we show the 
magnetoresistivity p^x obtained using our model, as a 
function of B for different MW field intensities, in all 
cases using the same frequency w/iir = v = 103.5 GHz. 
The darkness case is also presented. As the field inten- 
sity is lowered, the pxx response decreases to eventually 
reach the darkness response. In the inset it is possible 
to see the calculated pxx vs B~^, which is roughly peri- 
odic in B~^ in agreement with experiment. The minima 
positions as a function of B are indicated with arrows, 
corresponding to = j -I- ^. In the minima correspond- 
ing to j = 1, ZRS are found. Although the qualitative 
behavior of p^x as a function of B is very similar to the 
experimental one, the absolute value is smaller. It could 
be due to the smaller carrier density, or the simplified 
model for the electronic scattering with impurities that 
we have considered. Fig. 3 shows pxx versus B for differ- 
ent MW frequencies. The upper figure corresponds to a 
range of small values and the lower figure to large ones. 
In both cases ZRS are reproduced very clearly for j = 1. 
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FIG. 2: Calculated magnetoresistivity pxx as a function of 
B, for different MW intensities but for the same frequency 
V — 103.5GHz. The darkness case is also presented. In the 
inset we show p^x vs , which is roughly periodic in 
(T=1K). 
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FIG. 3: Magnetoresistivity response for different MW fre- 
quencies, a) corresponds to a range of low w values and b) to 
higher ones. In both cases, ZRS are reproduced (T=1K). 
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FIG. 4: pxx versus B for different temperatures with con- 
stant power excitation. The osciUations get smaUer as T is 
increased, but the positions of the minima stay constant. In 
the inset we show a calculated linear response between 7 and 
T. 

Minima positions shift as w is altered, maintaining the 
ratio — = 7 + T- This is shown in the inset, where the 
magnetic field B corresponding to ZRS (j=l) for four dif- 
ferent MW frequencies is plotted as a function of w. As 
in experimental results presented in ref. [1], a quite rea- 
sonable linear fit is achieved. 

The dependence of pxx on T (Fig. 4), has been ob- 
tained at = 103.5 GHz. As T increases pxx is soft- 
ened, and eventually almost disappears. The explanation 
can be readily obtained through the damping parameter 
7. When the electronic orbits are oscillating harmoni- 
cally due to the time-dependent external force, interac- 
tion with acoustic phonons occurs. This interaction acts 



as a damping for the orbits' movement. As T increases, 
the lattice-orbit interaction strengthens and the damp- 
ing of orbit dynamics will be stronger as well, giving a 
progressive reduction in the MW-induced pxx response. 
We have considered a linear dependence between 7 and 
T as in the experiments by Studenikin et al^. 
Using our model, it is now possible to shed some light 
on the peculiar dependence of pxx on B. According to 
our calculations we have found an approximately linear 
relation between pxx maxima and B, and we can there- 
fore express the corresponding dependence as: pxx 
Bcos{wt) (X i?cos (5). Looking at the cosine argument 
{wt) it is clear that if we change w, the minima posi- 
tions will change as well. Regarding r — — we can 
say that the scattering transition rate Wn,m is mostly 
dependent on sample and scattering variables, and that 
specific minima positions will be a function of those vari- 
ables. In this way we expect that for significantly dif- 
ferent samples, minima positions and other features of 
Pxx oscillations will change, explaining the discrepancy 
observed between different experiments. 
In summary, we have presented a new theoretical model 
whose main foundations are the exact solution for the 
quantum harmonic oscillator in the presence of MW radi- 
ation and elastic scattering due to randomly distributed 
charged impurities. This model gives a description of the 
electronic orbit dynamics which is crucial to explain the 
physical origin of ZRS. We are able to reproduce most 
experimental results, including pxx oscillations, minima 
positions and their dependence on w, MW intensity and 
T. 
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